Waveguide diffusion modes and slowdown of D'yakonov-Perel' spin relaxation in 

narrow 2-D semiconductor channels. 
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We have shown that in narrow 2D semiconductor channels the D'yakonov-Perel' spin relaxation 
rate is strongly reduced. This relaxation slowdown appears in special waveguide diffusion modes 
which determine the propagation of spin density in long channels. Experiments are suggested to 
detect the theoretically predicted effects. A possible application is a field effect transistor operated 
with injected spin current. 
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In connection to the possible use of electron spin for 
storage and information transfer in quantum comput- 
ers ||l| , there has been much recent studies on spin trans- 
port in semiconductor nanostructures. Among various 
materials for the so-called spintronic devices, a favorite 
candidate is III-V semiconductors because their spin- 
orbit split conduction bands have unusual electron spin 
dynamics. The spin-orbit interaction (SOI) in such ma- 
terials has the form Hso=h(k)-s, where s is the electron 
spin, and the direction and magnitude of the vector h(k) 
depend on the electron momentum k. When an elec- 
tron propagates, its spin precesses around the direction 
of h(k). In materials with narrow band gap, the mag- 
nitude |h(k)|, and hence the angle of spin rotation can 
be varied by simply applying a gate voltage 0. While 
this effect is valuable to those devices which make use 
ofballistically propagating spin-polarized electrons |g[ , it 
requires a high degree of coherency in electron propaga- 
tion, because scattering from an impurity or the bound- 
ary changes the direction of h(k). Scatterings thus ran- 
domize spin precession. If the mean scattering time r is 
shorter than the precession frequency |h(k)|~^, the spin 
dynamics is more like a diffusive relaxation in the an- 
gular space with the D'yakonov-Perel' relaxation rate of 
the order of l/rs=^r/i^(k) Q. However, the evolution 
of an inhomogeneous spin polarization at a spatial region 
is determined not only by the diffusive randomization of 
the local spin orientation, but also by a balance of the 
incoming and outgoung spin currents. These spin cur- 
rents depend on the polarization gradient, and contain a 
component associated to the spin precession. Hence, it 
is impossible to determine unambigously which part of 
the spin relaxation is due to the diffusion of spin orien- 
tation in the angular space, and which is due to the spin 
diffusion in the coordinate space. On the other hand, 
one can examine the evolution of a given spin polariza- 
tion by studying the dynamics of individual eigenmodes 
of the spin diffusion equation. After a sufficiently long 
time, only those modes with lowest spin relaxation rate 



survive, and this rate will be characteristic to the spin 
relaxation of the system under consideration. 

In this paper we will study the spin diffusion in a 2D 
channel made from a quantum well with growth direc- 
tion along the y axis. We choose our coordinate sys- 
tem to have x axis along the channel, and the bound- 
aries of the channel are marked at z=ztd/2. The channel 
width d is much shorter than the spin precession length 
Ls^Vp/lhCkp)], where Vp is the electron Fermi velocity. 
On the other hand, the width d is much longer than the 
Fermi wavelength, so that the electron motion perpen- 
dicular to the channel (along z axis) is semiclassical. We 
will show that the long-time evolution of spin polarization 
in such a channel is dominated by a waveguide diffusion 
mode. The spin relaxation time of this mode is slowed 
down dramatically by a factor L'^/d^ with respect to the 
D'yakonov-Perel' relaxation time Ts, which is the typi- 
cal spin relaxation time for bulk materials and quantum 
wells. Furthermore, the diffusion mode produces a peri- 
odic rotation of the spin polarization from point to point 
along the 2D channel. Such a phenomenon is similar 
to the spin precession of ballistic electrons considered in 
Ref. [Q , but with a difference that in our case the rotating 
spin has the spin quantum number 1 instead of 1/2. This 
difference is due to the fact that our nonequilibrium spin 
density is represented by two-particle excitations with 
electrons above the Fermi level and holes below it. As an 
example, we will demonstrate the oscillations in the re- 
sistance of a diffusive channel when the spin polarization 
is injected and probed via ferromagnets at both ends of 
the channel. We will also discuss qualitatively the effects 
of spin relaxation slowdown in a channel on the weak 
localization behavior of transport parameters. 

In terms of the creation and destruction operators cl^{t) 
and Ck(t), the spin density S(r,t) of a 2D electron gas is 
defined as 

S(r,t) = ^e''i--(Tr[c^^(i)sck(i)]), 

k,q 

where the average is taken over the ground state of 



the electron system. Taking into account the SOI and 
the less important Coulomb effects on spin density ex- 
citations, from the quasiclassical kinetic equation for 
<rr[c|^^ (f)sck(0]>, 01' by means of the standard per- 
turbation theory 0i^,^i one can derive the diffusion 
equation for the spin density S(r,t). 

It is convenient to represent S(r,i) in the basis set 
of three eigenstates of the z-component Jz of the an- 
gular momentum operator J which has the angular mo- 
mentum quantum number J—1. Accordingly, we intro- 
duce ?/'i=(5'a;-i5'y)/A/2, V'o='S'2, &nd'4)_i=-{Sx+iSy) / \/2, 
where the indices 1, 0, and -1 are the three eigenvalues 
of Jz- In term of this basis set, the diffusion equation is 
expressed as 



'di 



+ T{{-iv^ • Vr + hk • Jf)d^ri^ = /(r, t) , (1) 



where {...)dir is an angular average over the Fermi line 
and /(r, t) represents a possible source of spin oriented 
electrons inside the channel. The corresponding eigen- 
mode equation is simply 



r((-^v^ • Vr + hk • if)d^A = T^ ■ 



(2) 



The eigenvalue F is equal to the relaxation rate of the 
corresponding diffusion eigenmode. Eq. (|l|) is valid if 
V p\W r'>P\-^l / T and |/i(fci?)|^l/r. The second inequal- 
ity can not be satisfied for some high mobility InAs 
based quantum wells with a strong SOI, because in which 
v^/\h{kF)\<hm nm according to Ref. § and @. Eq. (0) 
can be generalized to the region \h{kp)\>\/T, a situation 
outside the scope of the present work. 

In a semiconductor quantum well, there are two con- 
tributions to the SOI h(k): the Dresselhaus term has 
its origin in the bulk crystal spin splitting p|, and the 
Rashba term is due to the asymmetric potential profile 
in the quantum well P]. The Rashba term is the domi- 
nating one in quantum wells made of narrow gap semi- 
conductors such as InAs [EJ. Therefore, for the systems 
of our interest, we will retain only the Rashba term of 
SOI, which can be expressed as hz=akx, hx—-aky, and 
hy—0. After taking the average over the direction of kp, 
Eq. (g) becomes 

D(i^- + m*aJzYiJ + Dii— ~ m*aJ^?^ = F?A , (3) 
ox oz 

where D^v^t/2 is the diffusion constant, and m* is the 
electron effective mass. If we replace D by l/2m*, the 
above equation is quite similar to the Schrodenger equa- 
tion for an electron. However, an important difference 
is that -0 is not a two-component spinor, but a three- 
component vector in the three-dimensional Hilbert space 
of the angular momentum J=\. 

^From Eq. (^, the spin flux can be expressed as 

F = -D{Vr + im*a3 x y)V' , 



where y is the unit vector along y axis. The boundary 
conditions at z=±d/2 are no spin flux across each bound- 
ary. Hence, 



d 
D{-i— + m*aJx)il^ \z=±i 



= 0. 



(4) 



For an extended 2D electron gas with d^oo, it is easy 
to see from Eq. (O) that the relaxation rates of homoge- 
neous spin density are given by the eigenvalues of the op- 
erator Dm*^a'^{Jz+Jx)- The two eigenvalues correspond 
to the transverse (spin polarized in xz plane) D'yakonov- 
Perel' spin relaxation rate l/Ts=Dm*'^a^, and the longi- 
tudinal (spin polarized in y direction) D'yakonov-Perel' 
spin relaxation rate 2/ts. In a 2D channel of electron 
gas, such homogeneous solution of Eq. (||) with spin po- 
larized along the channel, which satisfies the boundary 
conditions Eq. (^ also exists in the form of ■0±i=il/'\/2 
and ipo—0- This state relaxes with the rate 1/ts as in an 
extended 2D electron gas. 

However, we will prove in this paper that in a nar- 
row 2D channel there exist inhomogeneous diffusion 
modes with much lower spin relaxation rates. To derive 
these modes we will perform a canonical transformation 
'0=t/(z)^, where U{z)=exTp(-iJxz/ L^) with Lg=l/m*a. 
The operator U{z) transforms Eq. (||) to 



D 



-I- m aJz[z) 

ox 



1 2 



^j + D-^iP^Ti;, (5) 

oz^ 



with the boundary condition {d^p/dz)\z=±d/2=0i where 
Jz{z)—U^^{z)JzU{z). If d<^Ls, U{z) can be expanded 
in powers of the small parameter z/Ls- To the second or- 
der we get Jz(z)=Jz+Jyz/Ls-Jz{z/LsY /2 in which the 
last two terms will be treated perturbatively. The low- 
est order perturbation result gives the eigensolutions of 
Eq. (P) 



i'M,k,7n{x, z) = exp{ikx)Xm{z)'^ M : 



(6) 



where M and \1/m are eigensolutions of Jz with AI—±1 
and 0. The Xm{z) functions are X2^(z)=cos(27r2;n/d), 
and Xj^^j (z)=sin[7rz(2n-f l)/d] with integer n. 

The unperturbed relaxation rates are readily obtained 
from Eq. (||) as 



^ ALk,\ 



D{TTm/df + D{k - ML] 



-l^2 



(7) 



Hence, the modes with m^Q relax very fast. On the 
other hand, the modes with m=Q, which correspond to 
Xg(z)=constant, can have very low spin relaxation rates, 
which become zero for k=MLj^. However, corrections 
due to terms depending on dL^^ in Eq. (||) make these 
rates finite. Using the standard perturbation theory to 
the second order oi dL~^, we find 



^M,k,o — r^ 



A/,fc,0 



(2 - AP)d^ 

24r,L2 



(8) 



for \k-AILj^\<^Lj^ . The above equation indicates that 
although the relaxation rates of the modes with k=ML~^ 
become finite, their relaxation times retain much longer 
than the D'yakonov-Perel' relaxation time r^. 

There are three slowly relaxing diffusion modes with 
771=0, Af=l, 0, -1, and k=MLj^. The two modes for 
M = ±1 correspond to inhomogeneous spin distributions 
along the channel {x axis) with wavevectors k=^Lj^. 
Linear combinations of these modes give periodic spin 
density distributions S{x) with electron spins rotating 
in the plane perpendicular to z-axis. Hence, such a 
spin density has a finite projection along the channel. 
The third mode for M=0 corresponds to electron spins 
pointing along z-axis. The lowest relaxation rate is then 
achieved by an inhomogeneous spin density distribution 
along the channel. We notice that although for m =0 
the eigenfunctions given by Eq. (0) are independent of 
z, their corresponding spin distributions S(r) depend on 
z. This is because S(r) is expressed via the '0 function, 
which is obtained by applying the z-dependent unitary 
transformation U{z) to ipM,k,oix). 

We consider a stationary spin distribution in a chan- 
nel of finite length L^d, with an applied electric voltage 
V between xi—-L/2 at the left and Xr—L/2 at the right. 
This voltage can produce a magnetization in the channel, 
as was first proposed by Aronov ||l^] . The channel is con- 
nected with two ferromagnetic contacts, and a spin flux 
Fi is injected into the channel at x—xi. At x=Xr a spin 
flux Fr is then collected. The corresponding boundary 
conditions are 



d 



= F, 



Lr ■ 



(9) 



The stationary spin density in the channel can be de- 
rived from Eq. (||) with above boundary conditions. How- 
ever, our goal is to calculate the change of DC resistance 
associated with the spin transport, which can be mea- 
sured experimentally. To tackle this problem, we will 
use the formalism developed by Johnson and Silsbce [Q . 
The semiconductor-ferromagnet contacts are replaced by 
symmetric tunnel junctions with identical conductances 
G=Ae'^N{Ep)Vpt/2, where A is the channel cross section 
area, i<Cl the transmission probability, and N{Ef) the 
density of states at the Fermi energy. 

A DC current through the system will cause intcrfa- 
cial voltage drops AV/ across the left junction, and AV,. 
across the right one. Gradients AS'/ ,. of nonequilibrium 
magnetization across the tunnel junctions are also gen- 
erated. Following Ref. (ll|, the spin fluxes through the 
junctions can be written as 



Fi^r = -G [ {f]l,r/2e)AVl,r - {^/e^)ASl,r ] 



(10) 



where parameters rji and 7y,- with |?//.r|<l depend on the 
magnetizations of the ferromagnetic contacts, rji and 7y,. 
have same sign if the magnetic polarizations in contacts 



are parallel, and opposite signs if antiparallel. The first 
term in Eq (10) represents the spin injected by the ap- 
plied voltage, and the second term with ^~1 is the spin 
transport driven by gradients ASi^r- If in ferromagnets 
this magnetization relaxes sufficiently fast, ASi^r are de- 
termined mainly by the spin polarizations in the channel 
at xi and Xr, and so ASi=-ip{xi) and ASr—^'i^r)- Under 
stationary condition, the DC current is given by [nl| 



I ^ -G[AVLr - {VLr/2e)ASl,r] • 



(11) 



In this paper we will investigate the two cases that the 
magnetizations of the ferromagnet contacts are either po- 
larized along the z axis (perpendicular to the channel) or 
along the x axis (parallel to the channel). For each case, 
the magnetizations of the two contacts may be parallel 
or antiparallel. In terms of the momentum operator J, 
the magnetization is represented by the M—0 state if it 
is along the z axis, but by a linear combination of the 
Af=±l states if it is along the x axis. Hence, the corre- 
sponding injected fiuxes are represented by Fi^r=fi,r''Po 
and F/ ,.=//. r('0i-'0-i)7 respectively. Since L^d, the 
dominating contribution to the stationary spin distribu- 
tion comes from the modes with low relaxation rates. 
The corresponding exponents of the stationary solutions 
of Eq. (H) can be derived from the equation rjv/.fe.o=0. 
In this way, from Eq. (ph we obtain k=kM=MLj^±il^,j , 
where ^J^:/ = ^(2 - AP)/2M/Ll. In order to satisfy the 
boundary conditions Eq. (ph, we have to make a proper 
linear combination of the exp (x/Im) and exp (—x/Im) 
solutions. Furthermore, to demonstrate the spin dynam- 
ics along the channel, we average Eq. (H) over z. It is 
easily to see that the corrections which are linear in d/Ls 
vanish after the averaging. If we ignore the higher order 
corrections, we can represent "0 as a linear combination 
of the zero order eigenfunctions Eq. (0) with m=0 and 
k—kM- The spin density distribution is thus represented 
by 

V' = ^ 0M e*^^^/^-' [Am cos\i{x/Im) + Bm sinh(a;/ZM)] • 



M 



(12) 



iFromEq. (|) we derive Ao=(//-/.)L,/[2i?cosh(L/2/o)] 
and BQ—-{fi+fr)Ls/[2Dsmh.{L/2lQ)\ for the magnetiza- 
tion along z axis. Similarly, for the magnetization along 
X axis, the coefficients A±i and B±i are obtained as 

MLs ,^ MIL iML 

^^' = 777; — r^~^j' ^^p^^n—) ~ /'■ exp(-ir^— )j , 

2D cosh T=^ 2Ls 2Ls 



Bm = — 



MLs 



2D sinh 



AML iML, 
-jrUi exp(^^^— ) + /,. cxp( 



2£,, 



2L,, 



It is important to point out that although the bound- 
ary conditions require the spin polarization at x^±L/2 
to be along the x axis, within the channel the spin po- 
larization rotates in the xy plane due to the exponential 



factors exp (±zMa;/I/s) in Eq. (|T^. Since Ls<€.Im, the 
polarization makes many complete rotations within the 
spin relaxation length Ls- As it will be shown below, this 
can lead to an oscillation of the channel resistance. 

Substituting the expressions of Am and Bm into 
Eq. (I2h, one can find the spin densities at a;=±L/2, 
and hence the nonequilibrium magnetizations AS/.^. To- 
gether with Eqs. ( p^ ) and (^, we have a closed set of 
equations to solve numerically. However, it is important 
to illustrate in analytical form how the resistance depend 
on the spin transport through the 2D channel. For this 
purpose we consider a sample of length L:^Im, and hav- 
ing a low enough transmission probability t such that the 
parameter K=tlM/l^i, where I is the electron mean free 
path. From Eqs. ([lO|)-(|l^), we readily derive the spin- 
transport correction to the DC resistance as 



AR = {2K/G){T]f + Til) coth{L /Im) 

-{2K/G)riirir cos{M L / L s) / smh{L / 1 m) ■ 



(13) 



The factor 77/77^ in the second term at the right hand 
side is positive if the magnetizations of the two ferro- 
magnetic contacts are parallel, but negative if antiparal- 
lel. Furthermore, when the magnetization is along the 
channel, M=±l and so this second term oscillates as a 
function of L/Lg- The amplitude of oscillation decreases 
with increasing L/Im- However, due to the slowdown of 
the spin relaxation, we can manipulate the sample pa- 
rameters such that Im—L>Ls. Then, our theoretically 
predicted resistance oscillation can be observed experi- 
mentally with a gate to change the value of Lg- 

In our analysis above we have neglected the Dressel- 
haus contribution to the SOI, because in narrow gap sys- 
tems which we are interested in, the Rashba contribution 
dominates the SOI. In quantum wells the Dresselhaus 
contribution contains a linear term and a cubic term in 
electron momentum. It can be shown that the cubic term 
gives rise to an additional spin relaxation which is inde- 
pendent of the channel width d. Hence this term imposes 
a limit on the slowdown of spin relaxation. The corre- 
sponding relaxation rate is estimated to be insignificant 
for narrow gap quantum wells. However, in GaAs based 
quantum wells the cubic term is not negligible Q, and 
can wash out the slowdown of spin relaxation. 

Besides the classical spin diffusion, the spin relaxation 
slowdown also affects the weak localization corrections 
to transport parameters. In systems with strong enough 
SOI, the sign of such a correction to conductance is de- 
termined by the competition between the triplet and the 
singlet component of the Cooperon propagator IQ . The 
lifetime of the triplet is equal to the spin relaxation time. 
If the temperature is not very low, in a sufficiently narrow 



channel jlj] this time can increase and becomes compara- 
ble to the dephasing time of the singlet. Hence, the sign 
of the correction can change from positive to negative. 

The weak localization corrections to the spin diffusion 
coefficient and the D'yakonov-Perel' relaxation rate also 
contain contributions from the triplet and the singlet 
components of Cooperon [^ 14|. The spin diffusion coef- 
ficient is included in the factor k in Eq. (13). Therefore, 
by measuring the change of channel resistance AR with 
a controlled gate voltage or a weak magnetic field, the 
variations of spin diffusion coefficient can be investigated 
experimentally. 
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